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SIGNAL FOLDING COHERENT ACQUISITION FOR WEAK GLOBAL POSITIONING 

SYSTEM (GPS) C/A CODED SIGNAL 

RIGHTS OF THE GOVERNMENT 

The invention described herein may be manufactured and used by or for the 
Government of the United States for all governmental purposes without the payment 
of any royalty. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to global positioning system (GPS) 
satellite signal acquisition and more specifically to a faster computation algorithm for 
GPS P(Y) code and multiple blocks C/A code satellite signal acquisition. 

The nominal GPS operational constellation consists of 24 satellites that orbit 
the earth in 12 hours. The control segment consists of tracking stations located 
around the world. The GPS user segment consists of the GPS receivers and the user 
community. GPS satellite transmits specially coded satellite signals that can be 
processed by GPS receivers and provide the information to compute a user's velocity, 
time and position. 

The GPS satellites transmit two microwave carrier signals. FIG. 1 is a prior art 
drawing of GPS signals. The LI frequency (1575.42 MHz) shown at 100 carries the 
navigation message. The L2 frequency represented at 105 (1227.60 MHz) is used to 
measure the ionospheric delay by precise positioning service equipped receivers. 
Three binary codes shift the LI and/or L2 carrier phase. 

The Coarse Acquisition Code (C/A) shown at 102 modulates the LI carrier phase. 
The C/A code is a repeating 1 MHz Pseudo Random Noise (PRN) Code. This noise-like 
code modulates the LI carrier signal, "spreading" the spectrum over a 2 MHz 
bandwidth. The C/A code repeats every 1023 bits (one millisecond). There is a 
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different C/A code PRN for each GPS satellite. GPS satellites are often identified by 
their PRN number, the unique identifier for each pseudo-random- noise code. The 
C/A code that modulates the LI carrier is the basis for the civil standard positioning 
service (SPS). 

5 The P-Code (Precise) shown at 104 modulates both the LI and L2 carrier phases. 
The P-Code is a veiy long (seven days) 20 MHz PRN code. In the Anti-Spoofing (AS) 
mode of operation, the P-Code is encrypted into the Y-Code. The encrypted Y-Code 
requires a classified AS Module for each receiver. Therefore, only authorized users 
with cryptographic keys can use it. The P (Y)-Code is the basis for the precise 
JM> positioning service (PPS). 

The navigation message shown at 103 also modulates the LI -C/A code signal. The 
In Navigation Message is a 50 Hz signal consisting of data bits that describe the GPS 
m satellite orbits, clock corrections, and other system parameters. 

; ~ The C/A code and P(Y) code are code division multiple access (CDMA) systems 

;§5 where a pair of unique signals are assigned to each satellite in the GPS phase of the 
111 C/A code or the P(Y) code. The GPS receiver applies correlation to measure timing. 
M: The received signal is correlated with the locally generated replicas of the selected 
satellite's signal. This process is called acquisition. The traditional GPS receiver 
acquires this phase by continuous sliding, multiplication, and addition. This process 
20 is time consuming and is not conducive to miniaturized receivers. The C/A code is 
used in civilian GPS receivers and the military GPS receivers use both C/A code and 
P(Y) code. In general, the military receiver acquires the C/A code and transfers this 
timing to P(Y) code for tracking. However, if the military GPS receiver is under hostile 
environment and exposed to a strong jamming threat, the less vulnerable direct P(Y) 
25 acquisition becomes necessary. The present invention applies to both the C/A code 
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and the P(Y) code to reduce the calculations required to find a weak GPS signal 
through coherent acquisition. 

The conventional P(Y) code acquisition uses a time domain correlation 
approach as shown in Figure 2. For each satellite, this approach correlates 10ms of 
5 received sampled data (500,000 data points if sampled at 50 MHz), represented at 200 
with 200 locally generated replica, represented at 201. These replica are represented 

by 

r(m) = Pj{mAtlexpij2Kt k mMi (1) 

i> 

; ; p where At is sampling interval, Pj{mAt) is the sampled P(Y) code of satellite j ? m=0 ? 1 3 2, 
;fp 499,999 is a time index, and, fkis the center frequency of the locally generated 

ftj replica. To acquire the P{Y) code of the received signal from a targeted satellite, 200 
locally generated replica are correlated with 500,000 sampled points of the received 
ijjfS signal. If any of these 200 correlation result is above the threshold which is pre- 
f H determined by the correlation noise floor statistics, the code and the carrier frequency 
H : acquisition is completed, as is represented at 202. If none of the results is above the 
threshold, another 500,000 sampled data will be processed in the same manner, as 
represented at 203. This new 500,000 data set, represented at 204, only shifts one 
20 data point from the previous one. This process continues until either a signal is 
found or 2 ms of search range is exhausted. For ±lms of search range, the average 
amount of mathematical operations is 200x100,000 500000-point correlation, making 
the known approach a time consuming and energy consuming operation. 

To find a GPS C/A code signal, one need to find its carrier frequency and the 
25 beginning of the code. The C/A code is repetitive every millisecond, thus it is 
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adequate to search for 1 ms of data if the signal is reasonably strong. An acquisition 
method is required to search in both the time and frequency domains. For a high 
speed aircraft, the expected Doppler frequency range can be +/- 10 KHz. Therefore, 
the acquisition will search 20 KHz for the frequency. Let us assume that the input 
signal is digitized at 5 MHz. One millisecond of data contains of 5,000 data points. 

One way to perform acquisition is through circular correlation in frequency 
domain. In this approach the following steps are performed. 

(a) Generate a group of complex radio frequency (RF) signals according to 
the length of the data to be processed. This data length should be a multiple number 
of 1 ms. The frequency of this signal must be +/- 10 KHz around the carrier 
frequency with the frequency resolution equal to 1 over the data length . Perform 
steps (b) through (g) for each complex radio frequency signal generated in step (a). 

(b) Multiply the input signal by one of the complex RF signals. The 
resultant signal is very close to base band (or zero frequency). The operation can be 
expressed mathematically as 

si(n) = s(n) - rf(n) (Eq. 1) 
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where s(n) is the digitized input signal and rf(n) is the digitized complex radio 
frequency (RF) signal in time domain, which can be expressed as 



rf (n) = ei 2 *^ 



(Eq. 2) 



(c) Transfer s\ from time domain to frequency domain through the fast 
Fourier transform (FFT) to frequency domain as 
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Si(k) = FFT[si(n)] (Eq- 3) 

(d) Take the digitized C/A code of a certain satellite and make them to 
match the input data length, which is a multiple of ms. Call this signal c(n). 

(e) Take the FFT of c(n) and then take the complex conjugate as 
C(k) = {FFT[c(n)]}* 4 ) 



=M) Where * represents taking complex conjugate. 

% (f) Multiply C(k) and Si(k) and take inverse FFT. The result is 



tjj r(n) = IFFT[Si(k)C(k)] (Eq. 5) 

IJ5 where r(n) is the circular correlation result of Si(n) and c(n). 

]2 (g) Take the absolute value of r(n) and find the maximum. If the maximum 

is greater than the precalculated threshold, the value n of the maximum is the 
beginning of the C/A code and the RF frequency is the frequency of the carrier. 
20 Consider the description of circular correlation with 1 ms of data. One 

millisecond of data in time domain can generate 1 KHz resolution in the frequency 
domain. Therefore, the acquisition program will search the 20 KHz frequency range in 
1 KHz steps. There are 21 frequency bins to be searched from -10 KHz to +10KHz. 
This approach performs 5,000 point FFT and inverse FFT at each frequency. 
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Consider the description of circular correlation with 10 ms of data. Ten 
milliseconds of data digitized at 5 MHz generates 50,000 points. Also, ten 
milliseconds of data in time domain can generate 100 Hz resolution in the frequency 
domain. Therefore, the acquisition program will search the 20 KHz frequency range in 
5 100 Hz steps. There are 201 frequency bins to be searched from -lOKHz to + 10KHz. 
Following the discussion in the previous section, this approach performs 50,000 point 
FFT and inverse FFT at each frequency. This approach can find weaker GPS signals. 
One can see that this approach increases the complication of calculation 
tremendously. 

1Q For a GPS receiver with an outside antenna, the signal strength is sufficiently 

E strong such that 1 ms of data can be used for acquisition to find the signal. However, 
for weak signals, such as placing an antenna indoors, a long record of data is needed 

in for acquisition. The present invention can reduce calculations for acquiring GPS 
~~ signals using a coherent acquisition method. 

ft SUMMARY OF THE INVENTION 

fH When performing GPS signal acquisition, the strength of a signal received at a 

I,* GPS receiver having an outside antenna is generally sufficiently strong such that only 

1 ms of data needs to be used in acquisition to find the signal. However, for weak 

signals received at a GPS receiver, when the antenna is inside a building or vehicle, 
20 for example, a long record of data is needed for acquisition. The present invention 

reduces the calculation burden for GPS signal acquisition for a weak signal. 
It is therefore an object of the invention to provide a simplified coherent 

acquisition method for weak GPS C/A code. 
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These and other objects of the invention are described in the description, claims 
and accompanying drawings and are achieved by a simplified, weak GPS C/A code 
coherent acquisition method comprising the steps of: 

Receiving and digitizing a weak global positioning system C/A code data signal of 
a length longer than 1 ms; 

generating a complex radio frequency digitized signal as a local reference 

signal; 

first multiplying said digitized data signal from said receiving step with said 
digitized complex radio frequency signal; 

dividing an output of said multiplying step into N equal sections; 

adding data in N equal sections together from said dividing step; 

first applying a fast Fourier transform to an output of said adding step; 

considering 1ms of C/A code of a preselected GPS satellite; 

second applying a fast Fourier transform on the output of said considering 

step; 

taking a complex conjugate of an output of said second applying step; 
second multiplying an output of said taking step with output of said first 
applying step; and 

taking an inverse fast Fourier transform of the output of said second 
multiplying steps, an index of the maximum of said inverse fast Fourier transform 
being the initial phase of the C/A code. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a prior art description of GPS satellite signals. 

FIG. 2 shows a prior art conventional time domain approach for a GPS receiver. 
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FIG. 3 shows a frequency domain output of a fast Fourier transform of 10 
periodic C/A code. 

FIG. 4 shows a frequency domain output of a fast Fourier transform of 1 C/A 

code. 

FIG. 5 shows a block diagram representation of an input signal and 10 cycles of 
C/A code. 

FIG. 6 shows a frequency domain graph of an output from 1 ms of data at 3000 

Hz. 

FIG. 7 shows a frequency domain graph of output form 10 ms folded data at 
2800 Hz. 

FIG. 8 shows a frequency domain graph of output form 10 ms data at 2800 Hz. 

DETAILED DESCRIPTION 

The present invention simplifies the calculation of a coherent acquisition 
method for weak coarse /acquisition (C/A) code. The present invention may be 
described as a signal folding acquisition method because several 1 ms of data are 
folded and added. For purposes of simplicity, and in order to be comparably 
descriptive with conventional methods previously described, the signal folding 
acquisition method will be illustrated using 10 ms of data, however the present 
invention is operable with a wide range of signal record lengths. 

Using a ten millisecond record of data, the code function c(n) in Eq. 4 repeats 
itself 10 times. Recall, c(n) is the digitized C/A code of a certain satellite signal. The 
frequency domain result C(k) generated through the FFT of such a periodic function, 
c(n), is shown in FIG. 3. In FIG. 3 , the x-axis at 301 represents frequency bins and 
the y-axis at 300 represents power density. In FIG. 3, the output should be 50, 000 
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points but only 51 points (from 0 to 50) are shown. Only 51 points are shown 
because the C/A code repeats itself each 1ms and the total data length is 10ms. 

FIG. 4 illustrates signal acquisition employing the present invention using a 1 
ms record length of C/A code and the C(k) values are obtained from a 5,000 c(n) 
5 through FFT. In FIG. 4 the x-axis represents frequency bin, the y-axis represents 
power density and the symbol w * w represents the power density of the fast Fourier 
transform of the C/A code. In FIG. 4 six points indicated by * are shown. Comparing 
FIGs. 3 and 4, one can see that the partial result of the 50,000 point C(k), illustrated 
in FIG. 3, can be obtained from a 5,000 point C(k), illustrated in FIG. 4. In other 

iO words, properly padding the 5,000 point C(k) with zeros produces the same result as 

;g the 50,000 point C(k). 

[5 Since most the values in the 50,000 point C(k) are zeros, the corresponding 

ill terms in Si(k)C(k) are also zero. The acquisition method can be simplified as follows: 
f (a) Generate a complex RF signal of 10 ms long as shown in Eq. 2 with a 

i J5 frequency resolution of 100 Hz. To cover 20 KHz frequency range, 201 frequency bins 
lH are required. 

l2 (b) Multiply a GPS antenna received input signal by the complex RF signal 

from (a) and the result is 

si(n) = s(n) • rf(n) (Eq. 6) 

20 where s(n) has 50,000 points. 

(c ) Divide si into 10 equal sections and each has 5,000 points. Add all the ten 
sections together and the result is S2(n), which has 5,000 points. Take the FFT of 
S2(n), the result is 

25 S 2 (k) = FFT[s 2 (n)] (Eq. 7) 
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(d) Take the C/A code of a certain satellite of 1 ms long and call this signal 

ci(n). 

(e) Take the FFT of ci(n) and then take the complex conjugate as Ci(k) and 
then multiply with S 2 (k) and take the inverse FFT. The result is 

n(n) = IFFT[S 2 (k)Ci(k)] (Eq. 8) 

where ri(n) is the circular correlation result of S2(n) and Ci(n). In this approach, the 
frequency search range is still 201 bins. However, in each frequency bin the FFT and 
inverse FFT only perform 5,000 points, minimizing the calculation burden required in 
conventional approaches. 

The signal folding acquisition of the present invention may also be illustrated 
in the time domain. FIG. 5 shows 50,000 point input signal si(n), represented at 500, 
and 50,000 points (or 10 cycles) of C/A code in the time domain, represented at 501. 
Each section of the input signal Si(n) is correlated with the same C/A code. 
Therefore, mathematically, one can add all the 10 sections of the signal and correlate 
it with one cycle of C/A code. It should be noted that the input signal must be shifted 
to base band through the multiplication of rf(n). 

FIGs. 6, 7 and 8 illustrate the improved signal acquisition of the present 
invention in the frequency domain using varying records of data at varying 
frequencies. FIG. 6 shows the result of using 1 ms of data. In FIG. 6, the x-axis at 
601 represents time index and the y-axis at 600 represents power density. The 
frequency is at -3000 Hz because the frequency step is 1 KHz. The signal is not 
clearly recognizable. 
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FIG. 7 shows the result of using 10 ms of data. Similarly, in FIG. 7, the x-axis 
at 700 represents time index and the y-axis at 701 represents power density. In 
FIG. 7, the result is the first 5000 points obtained through 50,000 point FFT and 
inverse FFT. The frequency step is 100 Hz, the carrier frequency is at -2800 Hz. The 
signal at 702 is clearly shown. 

Finally, FIG. 8 shows the output using the signal folding acquisition method of 
the present invention. In FIG. 8 the x-axis at 80 1 represents time and the y-axis at 
800 represents power density. In FIG. 8, the result is obtained from 5,000 point 

FFT and it is almost identical to the result shown in FIG. 7. FIG. 8 illustrates that the 
folding method of the present invention requires less calculating but provides the 
same result as the conventional acquisition method for acquiring GPS C/A signals 
longer than 1 ms. 

While the apparatus and method herein described constitute a preferred 
embodiment of the invention, it is to be understood that the invention is not limited to 
this precise form of apparatus or method and that changes may be made therein 
without departing from the scope of the invention which is defined in the appended 
claims. 
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